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ABSTRACT 



We propose an explanation to the puzzling appearance of a wide blue ab- 
sorption wing in the He I A 10830 A P Cygni profile of the massive binary star 77 
Car several months before periastron passage. Our basic assumption is that the 
colliding winds region is responsible for the blue wing absorption. By fitting ob- 
servations, we find that the maximum outflow velocity of this absorbing material 
is ~ 2300 km s _1 . We also assume that the secondary star is toward the observer 
at periastron passage. With a toy- model we achieve two significant results. (1) 
We show that the semimajor axis orientation we use can account for the appear- 
ance and evolution of the wide blue wing under our basic assumption. (2) We 
predict that the Doppler shift (the edge of the absorption profile) will reach a 
maximum — 3 weeks before periastron passage, and not necessarily exactly at 
periastron passage or after periastron passage. 

Subject headings: (stars:) binaries: general— stars: mass loss— stars: winds, 
outflows— stars: individual (77 Car) 



1. INTRODUCTION 

There is no boring feature when it comes to the massive binary system rj Car. Its 
P = 5.54 yr (P = 2022.7 ±1.3 d; Damineli et al. 2008a) periodicity has been observed in 
the radio (Duncan & White 2003), IR (Whitelock et al. 2004), visible (e.g., van Genderen et 
al. 2006), X-ray (Corcoran 2005), and in many emission and absorption lines (e.g., Damineli 
et al. 2008a,b). Every orbital cycle the system experiences a spectroscopic event, defined 
by the fading, or even disappearance, of high-ionization emission lines (e.g., Damineli 1996; 
Damineli et al. 1998, 2000, 2008a,b; Zanella et al. 1984). The rapid changes in the contin- 
uum, lines, and in the X-ray properties (e.g., Martin et al. 2006, a,b; Davidson et al. 2005; 
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Nielsen et al. 2007; van Genderen et al. 2006; Damineli et al. 2008b; Corcoran 2005) are 
assumed to occur near the periastron passages of the highly eccentric, e ~ 0.9, binary orbit 
(e.g., Hillier et al. 2006). 

The He I A 10830 A high excitation line has a complex P Cygni profile, composed of 
three blue-shifted peaks with significant variations over the cycle (Damineli et al. 1998; 
2008b). The emission profile has significant variations over the cycle. The Doppler shifts of 
the peaks are of relatively low velocities, |w pea ks| < 300 km s _1 (Damineli et al. 2008b). The 
location of the minimum of the profile (the deepest point in absorption) does not change 
with orbital phase, and stays at v i >s _ m = —570 km s -1 . However, the absorption profile 
does change. In particular, just before periastron a wide blue wing appear in absorption, 
reaching ~ —1000 km s^ 1 a month before periastron, and ~ —1800 km s _1 at periastron, 
and the maximum equivalent width of absorption occurs 10 day after periastron passage 
(Damineli et al. 2008b). Damineli et al. (2008b) approximated the average radial velocity 
of the absorption profile at half intensity, and found it to change from —640 km s _1 before 
phase zero to —450 km s _1 shortly after phase zero. 

An additional He I A10830A blue absorption feature of up to —1000 km s -1 , is observed 
at several arcseconds from the center in the lobes (Smith 2002). In this paper we refer only 
to the He I A10830A ground measurements of Damineli et al. (2008b). We note that other 
absorption and emission lines of He I can be formed in different regions in the binary system 
(see also Kashi & Soker 2007b). In particular, some visible He I lines can be formed in the 
hot winds of the two stars close to their origin, as compared with the He I A10830A that 
is formed in cooler regions. We show that this cooler region can be the post-shock primary 
wind. Therefore, different He I lines need not have the same behavior along the orbit. 

Because of the winds' very complicated flow structure, when starting this project we 
limited ourself to build a toy- model in order to achieve two goals: (1) To show that the 
orientation where the secondary is toward us at periastron (u = 90°), can be accounted 
for the development of a wide blue absorption wing starting several weeks before periastron 
passage. (2) To encourage a nightly observation of the He I A10830A close to periastron 
passage. 

2. THE STARS AND THEIR WINDS 

The r) Car binary parameters used by us are compiled by using results from several 
different papers (e.g., Ishibashi et al. 1999; Damineli et al. 2000; Corcoran et al. 2001, 
2004; Hillier et al. 2001; Pittard k Corcoran 2002; Smith et al. 2004; Verner et al. 2005). 
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The assumed stellar masses are Mi = 12OM , M 2 = 30M Q , the eccentricity is e = 0.9, 
and the orbital period is P = 2024 day. The mass loss rates and terminal speeds are 
Mi = 3 x 1O" 4 M yr" 1 , M 2 = 1O" 5 M yr" 1 , v 1)0O = 500 km s~ l and v 2 ,oo = 3000 km s _1 . 
For these parameters, the half opening angle of the wind-collision cone is cj) a — 60° (Akashi 
et al. 2006); because of the orbital motion this angle is not constant, and we calculate it 
along the orbit in the present paper. For the inclination angle we take i = 41° (Davidson et 
al. 2001; Smith 2002, 2006). 

The primary's wind speed depends on latitude (Smith et al. 2003). The minimum 
in the He I A10830A line profile suggests that the primary's wind speed toward us is 
tUs-m = —570 km s _1 . The two winds collide, and form a flow structure, schematically 
drawn in Figure [TJ The two winds go through two respective shock waves, and form a 
contact discontinuity between them. The contact discontinuity asymptotically forms a con- 
ical shell surface. The radiative cooling time of the post-shocked primary's wind is short. 
The post-shocked primary's wind forms a dense flow along the contact discontinuity, which 
we refer to as the conical shell. Absorption is expected to take place mainly outside the 
stagnation-point region, and so we approximate the conical shell as an ideal cone. 

The He I A10830A absorption profile results from He I atoms which absorb from the 
He I A10830A emission line and from the continuum emitted by dust. As we are more 
interested in the wide wing, the continuum is more relevant to our study. As the system 
approaches periastron, dust is formed closer and closer to the binary system (Kashi & Soker 
2008a). This makes things very complicated, as the conical shell, where the absorbing gas 
reside according to our model, also changes its size. Some of the dust is formed in the conical 
shell itself. 

Other lines, e.g. Ha (Smith et al. 2003), also show some fast blueshifted absorption 
wings. Some of these wings disappear when the system approaches periastron. This is 
not a problem for our model, since those lines originate in different regions than where the 
He I A10830A line does (e.g. the Ha comes from the primary's wind; Smith et al. 2003). 
We emphasize again that some He I lines (in particular in the visible band) are formed in 
different regions than the He I A10830A line, and their behavior is different as well. 

The binary parameter that is most controversial is the orientation of the semimajor 
axis— the periastron longitude. Some researchers argue that the secondary (less massive) 
star is away from us during periastron passages, i.e. an orbital longitude of u — 270° (e.g., 
Nielsen et al. 2007; Damineli et al. 2008b), others argue that the secondary is toward us 
during periastron passages, lu = 90° (Falceta-Gongalves et al. 2005; Abraham et al. 2005; 
Kashi & Soker 2007b, 2008b), and still different values exist in the literature (Davidson 
1997; Smith et al. 2004; Dorland 2007; Henley et al. 2008; Okazaki et al. 2008a). Following 
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Fig. 1. — The colliding winds structure. The primary's and secondary's winds mix close to the contact 
discontinuity. According to our models this region is the absorber of the wide He I A10830A blue wing. The 
half opening angle of the wind-collision cone is <f) a ~ 60°. The thick arrows indicate the flow lines. 

our recent paper (Kashi & Soker 2008b) we will take the orientation to be such that the 
secondary is toward us at periastron {u = 90°). 

3. THE TOY-MODEL 

We suggest that the conical shell is responsible for the absorption of the He I A10830A 
high excitation line. When the two winds meet at the contact discontinuity, part of the 
shocked fast secondary's wind is mixed with part of the shocked slow wind (Pittard 2007), 
and accelerates part of it to higher velocities. The conical shell has an asymptotic angle of 
<f) a ~ 60°, which we take to be the flow direction of the absorbing gas. 

In order to make the model more accurate we will take into account the time dependance 
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of some parameters. The primary's wind velocity profile can be described using the /3-profile: 

/ R \ ^ 

vi(r) = v s + (vi t0O - v s ) ( 1 -) , (1) 



r 

where v s = 20 km s -1 is the sound velocity on the primary's surface, Vi j00 = 500 km s _1 is 
the primary's wind terminal velocity, and (5 = 1 is a parameter of the wind model. 

The radial (along the line joining the two stars) component of the relative velocity 
between the secondary star and the primary's wind is v% —v r , where v r the radial component 
of the orbital velocity; v r is negative when the two stars approach each other. The total 
relative speed between the secondary and the primary's wind is 

1 /2 

^windl = [Vq + Ol - V r ) 2 ] , (2) 

where vq is the tangential component of the orbital velocity. The orbital motion and the 
variation of the primary wind speed with distance from the primary have a small influence 
on the conical shell asymptotic angle <fi a . We will use the expression given by Eichler & Usov 
(1993) 




a ~2.l|l--V)/A (3) 
where 

, s [a^Z. (4) 

V Miw wind i 

We will take into consideration the rotation of the cone relatively to the line connecting 
the two stars. This rotation occurs due to the orbital velocity of the conical shell, and has a 
considerable influence close to periastron. We define 5(f) to be the angle measured from the 
secondary between the direction to the primary and that to the stagnation point (see Soker 
2005 for further details) 

cos ( 5 0) = V 1^JL ( 5 ) 

^windl 

We find that close to periastron 5<p — 56°. The geometry and different parameters are shown 
in Figures \T\ [2] and El 

A closely-related geometry was used by Hill et al. (2000) to fit the excess emission 
observed in the C III A5696A line in the spectra of WR 42 and WR 79. Luehrs (1997) also 
used a somewhat different geometry to fit the excess emission observed in the C III A5696A 
line in the spectra of WR 79. These two models were purely geometric. In the model we 
suggest for the absorption, some more considerations have to be taken in to account. 
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Fig. 2. — The geometry of the toy-model and definition of several parameters. The geometry is plotted 
near periastron (9 = 0), where according to our model the secondary is toward us. An example of a "tube" 
is also shown. Both 9 and 6<f> are measured in the equatorial plane. At each point along the orbit there is 
one value for each of the variables 9, 5<f>, and </) a , while we integrates over many tubes on the conical shell, 
each with its value of direction to the observed if). 

Although the absorbing gas is a continuous media, we decompose the conical shell into 
'tubes'. We define ( to be the azimuthal angle on the surface of the conical shell; ( = in 
a direction perpendicular to the equatorial plane, and is measured clockwise. Tubes exist 
from ( = to ( = 2tt. For each tube on the conical shell we calculate the angle if) between 
the tube and the observer as a function of the orbital angle 9 (see Figure [2D 

cos if> = cos 4> a sin % cos(# — 8(f)) + sin (f> a sin ( sin i sin(# — 5(f)) + sin (f> a cos ( cos i. (6) 

We take v o < v m < v max to be the maximum velocity attended by a mass element, 
dm, in a tube; in each tube there is a large number of mass elements formed by the mixing 
process of the two winds. By mixing we refer also to secondary shocked wind segments 
that have been cooled to temperature low enough to contribute to the absorption at high 
speeds. They reside near the contact discontinuity as well. Here vq = Vi >00 = 500 km s" 1 is 
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Fig. 3. — The time variation of angles defined in the text. 9 is the orbital angle (line from primary to 
secondary) with 9 = at periastron. 5(f) is the angle between two lines from the secondary: one to the 
primary star, and one to the stagnation point of the colliding winds. Both 9 and 5(f> are measured in the 
equatorial plane. <f> a is the opening angle of the cone, defined around the cone axis. 

the primary's wind velocity, and v max is a parameter of the model that is constraint to be 
i; max < t>2,oo = 3000 km s _1 . The projected velocity of mass element having a velocity v m 
and residing in a tube with an angle of ip to our line of sight is vd = v m cosi/j. 

Consider one of our tubes; it has a length L and a circular cross section with a radius 
Rt, such that Rt <C L. The effective cross section of the absorbing tube is 2R t Lsinifj. 
Therefore, the contribution of each tube to absorption is multiplied by sinip. The amount 
of accelerated primary's wind to each velocity v m is hard to predict. A numerical simulation 
beyond the scope of this work needs to be done in order to determine this amount. We 
simply take the velocity distribution to be constant, namely, in each tube the fraction of 
the gas in the velocity interval dv is constant W(v)dv = Cwdv m , where Cw is a constant. 
Changing variable from the velocity along the tube to that along the line of sight, this 
weighting function (up to a constant) is, by using vd = v m cosip, 

C 

W(v m )dv = W(v D )dv D = -dv D , (v cosip) < v D < (v max cosip). (7) 



-8- 



At each orbital angle 9 we sum the contribution to absorption at Doppler shift vp over all 
tubes, 



where K a is a constant of the toy-model that takes care of units and the condition < / < 1. 
We set herein = l/max(A), because we do not compare to the absorption equivalent width. 

Our model is actually a toy-model. It assumes a simple geometry of the absorbing 
material, i.e., a rotated conical shell with a varying opening angle, and a mass distribution 
within each tube that is constant with the velocity. The model contains two types of param- 
eters: (1) Those that are given in the literature. Such are the binary parameters and the 
conical shape with its opening angle (f) a . These parameters are more or less in consensus. 
The orientation of the semimajor axis (the periastron longitude u) is controversial, and we 
take the value from our previous papers u> = 90° (secondary toward us at periastron). (2) 
Parameters that are unique to our toy-model. Such is the value of t> max which is constraint 
to be v max < t>2,oo = 3000 km s _1 . For these parameters we find that a good general fit can 
be obtained for f max = 2300 km s _1 . The value of the absorption coefficient K a has a small 
influence on our conclusions, and it serves only to give the general form of the absorption 
profile. 

The assumption that the conical shell reaches its asymptotic opening angle at large 
distances breaks down as the system approaches periastron. The reason is that the relative 
velocity of the two stars is no longer much smaller than the primary's wind speed. This 
causes the winding-up of the conical shell into a spiral structure in the equatorial plane, e.g., 
as shown in the numerical study by Okazaki et al. (2008b). Only a close region near the 
binary system reaches this limiting angle. Namely, a smaller region, but much denser, will 
contribute to the absorption at the blue edge. However, contribution to the continuum near 
1 /im comes from the stellar wind and hot dust from closer regions to the binary system 
(Kashi & Soker 2008a). Therefore, the much smaller region of the conical shell can still 
absorb a detectable fraction of the continuum. Our toy-model does not allow us to make 
any quantitative prediction. For that a 3D numerical code is required. Nevertheless, our 
toy-model does allow us to reach the two goals, as mentioned in the last paragraph of section 





tubes 



The intensity is taken to be 




(9) 



X 
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4. THE BLUE ABSORPTION WING 

In Figure H] we present our results by a contour map of the intensity /, as given in 
equation EH in the Doppler shift— time plane. The blue Doppler shift Vd is given in unit 
of km s _1 , and the vertical axis indicates days relative to (before) periastron (phase at 
t = 0). The levels of the contours are I = 1,0.95,0.9,0.85,0.8,0.7,0.6, and so on, from 
left to right, where 1 = 1 indicates the edge of the absorption profile. Despite the constant 
weighting function W that we used, in reality we expect that the amount of helium blown 
with high velocities corresponding to I ~ 0.9 — 1 will be too small to be detected. We have 
also estimated the velocity edges of the absorption wing at four epochs from Damineli et al. 
(2008b) observations; these are indicated by four horizontal error-bar lines,. The somewhat 
noisy data made it difficult to pinpoint the exact edge of the wing, and therefore we could 
only estimate it to an accuracy of ~ 100 km s^ 1 . Each line is centered at the approximated 
edges of the absorption wing, and extends to ±50 km s" 1 to each direction. At early times 
there is no noticeable differences between the contour lines in the range / = 0.85 — 1, and all 
nicely fit the two observation equally well. Very close to periastron, the line / = 0.85 fits the 
observations better. However, during this time we expect the collapse of the conical shell, 
such that no fresh gas will be accelerated to high velocities. We cannot make an accurate 
prediction so close to periastron passage. Over all, we consider the lines / = 0.85 — 1 to be 
a very good fit to observations. 

We modified the parameters in our model to fit the observations of Damineli et al. 
(2008b) from the 2003 event. With this fitting we learn about two properties of the models. 
(1) We find that the orientation we use, where the secondary is toward us at periastron 
(u = 90°) can account for the development of a wide blue absorption wing starting several 
months before periastron passage, if the absorbing material is within this conical shell. (2) 
According to our model, the maximum observed blue shift of the wing is reached ~ 5 day 
before periastron passage. The observed maximum blue shift might occur at a somewhat 
different time for three reasons. Firstly, the binary and winds parameters might be somewhat 
different than those we used. Secondly, we did not take into account the winding of the conical 
shell; we only considered the approximate direction of its axis by calculating the angle 5(f). 
Thirdly, close to periastron the conical shell is likely to collapse onto the secondary (Soker 
2005; Kashi & Soker 2009). 

With the only four available observations for the fitting, with the present possible de- 
gree of accuracy, and according to our parameters, we expect the maximum blue shifted 
absorption to occur — 20 day before periastron passage, i.e., late December or early Jan- 
uary. There are two main reasons for that we cannot be more accurate. Firstly, we cannot 
treat properly the conical shell as the system approach periastron. Secondly, we expect the 
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Fig. 4. — A contour map of the intensity I as given in equation [51 in the Doppler shift— time plane. The 
blue Doppler shift is given in unit of km s , and the vertical axis indicates days before periastron (phase 
0). The levels of the contours are 1 = 1, 0.95, 0.9, 0.85, 0.8, 0.7, 0.6, and so on, from left to right. The edges 
of the blue absorption wing at four epochs from Damineli et al. (2008b) are marked. The length of each of 
the drawn lines corresponds to 100 km s^ 1 , our estimate of the uncertainty in determining the edge of the 
absorption. We predict the edge Doppler shift to reach a maxima 0-20 days (according to our parameters) 
before periastron, when a maximum in t)£> appears according to our model. 

colliding wind region to collapse onto the secondary near periastron; we did not consider this 
process here either. 

Our fundamental assumption is that the absorber of the blue wing of the He I A10830A 
line resides in the conical shell formed by the colliding winds. This assumption works quite 
well with the semimajor orientation u = 90°, where the secondary is closest to the observer 
at periastron passage. We now show that other orientations that are popular in the literature 
cannot account, not even qualitatively, for the behavior of the blue wing, if our fundamental 
assumption holds. We followed the calculations presented in Figure HI but for three different 
semimajor axis orientations, as drawn in Figure [5j uo = 0: The semimajor axis is perpendicu- 
lar to the line of sight and the secondary is closer to the observer before the event, to = 180°: 
The semimajor axis is perpendicular to the line of sight and the secondary is closer to the 
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observer after the event, uo = 270°: The secondary is closest to the observer at apastron, 
opposite to our favorite orientation of uj = 90° (Figure HJ). 
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Fig. 5. — The fours semimajor orientations studied in the paper. In all cases the inclination angle (the 
angle between line of sight and a line perpendicular to the orbital plane) is i = 41°. Our favorite model has 
u) = 90°, for which the blue absorption wing is drawn in figure [?J 

The results for the expected Doppler shifts of the absorption wing are presented in 
Figure El together with the fours observations from Damineli et al. (2008b). It is clear 
that none of the other orientations can account for the observed absorption wing, not even 
qualitatively and just for the two early observations. 



5. SUMMARY AND PREDICTION 

We study the blue absorption wing of the He I A10830A P Cygni profile of rj Car. The 
two winds of the two stars collide, and the post shocked gas of the two winds flow on both 
sides of a surface (the discontinuity surface) that has a pseudo-conical structure: the conical 
shell. The shocked secondary's wind accelerates part of the shocked primary's wind to high 
velocity (Figure [T]). This gas, and the segments of the shocked secondary's wind which cool 
to a low temperature near the contact discontinuity (Figure H]), are assumed to be responsible 
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for the blue absorption wing. This is the fundamental assumption of our model. We use our 
previous results and assume an orientation where the secondary is toward the observer at 
periastron (u = 90°; see Figure [5]). For the conical shell we built a toy model (Figure [2]). 

The absorption profile, up to a scaling factor, is calculated according to equations (jHj) 
and ©, and the results are presented in Figure HI We are interested in the bluest part of 
the absorption profile, where intensity is lower; in our scaled units these are the contours in 
the range / ~ 0.8 — 1. Using our fundamental assumption, and a toy model for the conical 
shell of the colliding winds, we showed that with our orbital orientation we can account for 
the appearance of the wide blue wing several months before periastron. Other semimajor 
orientations u, cannot reproduce the results under our fundamental assumption (figure E]) 
This is our main result, namely, that if the absorber responsible for the blue wing of the 
He I A10830A line reside in the winds collision region, then only the u ~ 90° can account 
for the blue wing of this line. 

Our results also predict that the Doppler shift vp (the edge of the profile) will reach a 
maximum — 3 weeks before periastron passage. Since close to periastron the conical shell 
starts to collapse onto the primary, and even before it experiences the effect of wrapping, it 
is hard to pinpoint the exact time of this maximum. Nevertheless this maximum should be 
observed before the event. 

We thank Augusto Damineli for providing us the observational data for the He I A10830A 
line, and an anonymous referee for helpful comments. This research was supported by the 
Asher Space Research Institute in the Technion. 
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Fig. 6. — Same as in FigureHl but for different values of w, as schematically depicted in Figured]. None of 
the other orientations can account for the observed absorption wing under our fundamental assumption that 
the absorber reside in the conical region formed by the collision of the two winds. In particular, the opposite 
orientation, lu = 270°, yields an entirely opposite behavior to observations, with minimum absorption at 
periastron. 



